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Abstract 

We study a holographic cosmological model in which the infrared cutoff is set by the Ricci's length 
and dark matter and dark energy do not evolve on their own but interact non-gravitationally with 
each other. This greatly alleviates the cosmic coincidence problem because the ratio between both 
components does not vanish at any time. We constrain the model with observational data from 
supernovae, cosmic background radiation, baryon acoustic oscillations, gas mass fraction in galaxy 
clusters, the history of the Hubble function, and the growth function. 
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I. INTRODUCTION 



The efforts made to unveil the nature of dark energy (DE) -the mysterious agent behind 
the present phase of cosmic accelerated expansion- have been scarcely rewarded so far. We 
only know by certain that DE is endowed with a hugely negative pressure (of the order of its 
energy density) and we have grounds to suspect that it is distributed rather evenly across 
space -see [l| to learn about the state of the art. 

In view of its mysterious nature many authors have suggested that dark energy should 
comply with the holographic principle which asserts that the number of relevant degrees of 
freedom of a system dominated by gravity must vary as the area of the surface bounding 
the system Qj. In addition, the energy density of any given region should be bounded by 



that ascribed to a Schwarzschild black hole that fills the same volume jj}]. Mathematically 
this condition reads px < MpL~ 2 , where px and L stand for the DE density and the size 
of the region (or infrared cutoff), respectively, and Mp = (87rG) _1//2 is the reduced Planck 
mass. This expression is most frequently written in its saturated form 

3M 2 c 2 

PX = -^2~ ■ (1) 

Here c 2 is a dimensionless parameter -very often, assumed constant- that summarizes the 
uncertainties of the theory (such as the number particle species and so on), and the factor 3 
was introduced for mathematical convenience. An interesting feature of holography lies in its 
close connection to the spacetime foam 4J. For additional motivations of holographic dark 
energy see section 3 of [5(]. At any rate, it is sobering to bear in mind that the holographic 
proposal is just a reasonable hypothesis (which we adopt in this paper) but not necessarily 
a compelling one. 

When dealing with holographic DE one must first specify the infrared cutoff. In the lack 
of a clear guidance different expressions have been adopted. The most relevant ones are the 
Hubble radius, i.e., L = H -1 , see e.g. {(J- and the Ricci's length, i.e., L = (H + 2H 2 )~ 1 / 2 



-see e.g. 



— Il0j| - The rationale behind the latter is that it corresponds to the size of the 



maximal perturbation leading to the formation of a black hole . The radius of the future 
event horizon have been profusely used but it suffers from a severe circularity problem. 

The aim of this paper is to present a cosmological model of holographic dark energy 
and constrain it with observational data. The model takes the Ricci's length as infrared 



cutoff and assumes that dark matter (DM) and dark energy do not evolve separately but 
interact non-gravitationally with one another. The interaction, albeit proposed just at a 
phenomenological level, is key to alleviate the coincidence problem 12|, [13J which cannot 



find explanation in the ACDM model and afflicts so many models of evolving DE. On the 
other hand, recently it has been suggested that the dynamics of galaxy clusters can be 
better explained if the said interaction is taken into account , and it has been invoked to 
explain the nearly lack of correlation between the orientations of cluster galaxy distributions 
with those of the underlying dark matter distributions 1^|. Further, the interaction -first 



proposed to lower down the theoretical value of the cosmological constant 
necessary but inevitable [17 ]. 



16] - is not only 



This paper is organized as follows. Section [TT] introduces the model. Section |TTT] discusses 
why Ricci's holographic models seems to solve the coincidence problem also in the absence 
of interaction. Section HVl presents the statistical analysis of the model after constraining it 
with data from supernovae type la (SN la), the shift of the first acoustic peak of the cosmic 
background radiation (CMB-shift), baryon acoustic oscillations (BAO), gas mass fraction in 
galaxy clusters, the history of the Hubble parameter, H(z), and the growth function. Lastly, 
section |V] briefly delivers our main conclusions and offers some final remarks. As usual, a 
zero subscript means the present value of the corresponding quantity. 



II. THE HOLOGRAPHIC MODEL 

This model assumes a spatially flat homogeneous and isotropic universe dominated by DM 
and DE (subscripts M and X, respectively), the latter obeying the holographic relationship 
([1]). In virtue of Friedmann equations, 

H 2 = l -Mf{p M + p x ) , (2) 

and 

H = - l -M p \p M + px+Px), (3) 
the fractional DE density, fix = ^1^772, can be expressed as 



where w = px/px denotes the equation of state parameter of dark energy which in general 
depends on time. 

The deceleration parameter q = —d/(a 2 H 2 ) takes the simple expression, 

(6) 

As we shall see, as a consequence of the evolution of fix, it goes monotonously from positive 
values at early times (in the matter dominated era) to negative values at later times (in the 
dark energy dominated era). 

The evolution of fix follows from the conservation equations of DE and DM. In the 
absence of non-gravitational interactions between them they evolve independently and obey 

nM-U-^r){i-n x )H = o, (6) 
n x +(i-^)(i~ nx)H = - (7) 

Bearing in mind that in our case + fix = 1, we get the following expressions in terms 
of the redshift (z = a" 1 — 1), 

n 2^o-c 2 + c 2 (l-Qxo)(l + #- 1 
2Qxo -c 2 + 2(1 - Q X0 )(l + 

A very useful quantity when considering the coincidence problem is the ratio between 
the energy densities, r = Pm/px, which for spatially flat universes reduces to fl M /fl X - The 
coincidence problem gets alleviated if for reasonable values of flxo and c 2 we get r ~ Oil). 
This is the case here since for flxo an d c 2 of order unity, r also results of this same order. 
However, for late times (i.e, when z — > —1) one has r — > 0. In this sense, the coincidence 
problem is not properly solved. To obtain a non- vanishing ratio at late times some interaction 



between DM and DE must be incorporated in the picture [18] . 

If DM and DE interact non-gravitationally with each other the evolution equations may 
be generalized as 

^M- (i-^r) (l-n x )H = QH, (9) 
flx+ (i-^r) (l-n x )H=-QH, (10) 

where the Hubble factor on the right hand sides has been introduced to render the interaction 
term, Q, dimensionless. 



Since the nature of both dark components is largely unknown, there is ample latitude in 
choosing Q. We shall specify it by demanding that r evolves from an unstable fixed point in 
the far past, = r(z — > oo), to a stable fixed point at the far future, 77 = r(z — > —1) [19 ]. 

The pair of equations (jU]) and ( ITUj) imply 



r ( 1 + r- % 1 + Q(l + r) 



c 2 



H. (11) 



Imposing that 77 be a fixed point, i.e., r\ r=rf = the interaction term Q is simply a constant 
given by 

« = -(rM 1+ "-?)- (12) 

As we shall see later, 77 and c 2 take values such that Q is positive-definite, which entails 
a transfer of energy from dark energy to dark matter. Obviously if Q were negative, the 
transfer of energy would go in the opposite direction which would conflict with the second 
law of thermodynamics [20] and the coincidence problem would only worsen. 
Rewriting Eq. ( ITT]) as 

r = (Q + l)(r-r f )(r-r oa )H (13) 

and using the condition r = 0, the other fixed point can be expressed in terms of the previous 
one, namely, 

- - 2 " c2 ( 1 + ^) (14) 



2r / + c 2 (l + r / )' 

To study the stability of the fixed points we first write r' = dr/dlna = r /H and calculate 
the derivative of r' with respect to r. In the case of the far future fixed point we get 

dr' _ 2(77-1) 

dr lr '~ 1+ 02(77 + 1)" (15) 
Since 77 must be lower than ro — 0.45, from Eq. (J5]) with c 2 < Qxo — 0.75 (otherwise qo 
would not be negative), one follows that ^r\ rf < 0, i.e., the fixed point 77 is a stable one. 
Similarly, we find that 

dr'. 2 + c 2 4 

^ = -— + ^T7J) > ^ (16) 

i.e., the fixed point at the far past is an unstable one. 

Equation fTT3|) can be integrated with the help of (TT2~]) . In terms of the redshift it yields, 

_ 77(7-0 - r^) - r^ro - r/)(l + z )(Q+ 1 )('--'7) 



(r - t-oo) - (r - r/)(l + z)^ +1 )^" r s) 



(17) 
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Inspection of ( Tl3|) readily shows that when r lies between both fixed points one has r < 0, 
i.e., the ratio between the energy densities diminishes monotonously from one fixed point to 
the other. This is depicted in Fig. [TJ The said ratio smoothly decreases from high z (i.e., 
from -the unstable fixed point is at z — > oo)) to asymptotically approach the fixed stable 
point, rf, at z — —1. Note that the latter needs not be zero. In this regard the coincidence 
problem is much alleviated because we are not living in any special era. However, the 
problem is not solved in full since the model cannot predict that tq is of order unity. To 
the best of our knowledge, no model predicts that, as well as no model predicts the present 
value of the temperature of the cosmic background radiation, the Hubble constant, or the 
age of the Universe. For the time being, we must content ourselves by taking these values 
as input parameters since, very possibly, we are to wait for a successful theory of quantum 
gravity to compute them. 

The expression for the fractional density of dark energy follows from the relationship 
r = (1 - n x )/ttx and Eq. (HTJ), 

(rp - roD ) - (To - r / )(l + gjW+lXrco-tT) 
X (r f + l)(r -r 00 )-(r 00 + l)(ro-r f )(l + z)(Q+ 1 ^- r f) ' 1 ' 

From the latter and (jl]) we obtain the equation of state of dark energy in terms of the 
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Figure 1. The ratio r between the energy densities vs. redshift for the best fit model. As the inset 
shows rf = r{z — > —1) does not vanish. In this, as well as in subsequent figures, the red swath 
indicates the region obtained by including the la uncertainties of the constrained parameters used 
in the calculation. 
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redshift, 



w 



1 / 2 

- [ 1 - - + ^ + 



{r - ^oo)(^oo - r f ) 



roo-r + (r - r/)(l + ;s )«+ 1 )(»-«>-'v) 



(19) 



As shown in the left panel of Fig. [2] for the best fit model, w smoothly evolves from a 
negative value close to zero at high redshifts to a value lower than —1 at the far future. The 
right panel depicts its evolution near the present time (z = 0) showing compatibility with 
recent observational data which suggest that w does not depart much from —1 at sufficiently 
low redshifts (see 12 IIP. 
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Figure 2. The equation of state parameter (as given by Eq. ([19]) ) vs. redshift up to z = 8 (left 
panel), and up to z = 1.2 only (right panel) for the best fit holographic model. At high redshifts w 



approaches the equation of state of non-relativistic matter and at 



significantly from —1. The observational data are taken from 
uncertainty. 



ow redshifts it does not depart 



211 ] : each error bar signifies a 2a 



Integration of the second Friedmann's equation ([3]) provides us with the evolution equa- 
tion for the Hubble factor which is key to perform the statistical analysis of section IIV[ 



where 



H = Hn 



A 1 + 2(A 2 + (r f -r )(l + z)- A *) 



A, 



1/2 



(l + zy 



A 1 = c 2 (l + r / )(l + r ), A 2 = r r f -1, A 3 = 1 + 
A 4 = [c 2 (l + Tf ) + 2(77 - 1)](1 + r ) , A 5 = 2 - 



2(r/ 



c\l + r f ) 
1 



(20) 



c=»(l + r/)" 

Figure [3] depicts the evolution of the deceleration parameter, Eq. ([5]), for the best fit 



model. The observational data are borrowed from 



22] . The redshift at which the universe 



la 
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Figure 3. The deceleration parameter vs. redshift for the best fit holographic model (solid line) and 
the ACDM model (dashed line). In this, as well as in subsequent figures, the green swath indicates 
the region obtained by including the la uncertainties of the constrained parameters used in the 
calculation (in the present case just £Imo)- The observational data are taken from 
in view of the large error bars we do not use these data to fit the models. 



22j |; however, 



starts accelerating is z(q = 0) = O.56T0/9 while for the ACDM model z(q = 0) = 0.79 ±0.02. 

As mentioned above, there is ample freedom in the choice of the interaction term Q. In 
a previous paper of two of us, on a holographic dark energy model with the Hubble rate as 
infrared cutoff, we took Q oc fix [23(. We do not pursue this possibility here because, as we 
have checked, it leads to a universe in which dark energy is subdominant at very late times. 
While this does not contradict observation, it looks a bit odd. In any case, it deserves a 
separate study which lies beyond the scope of this paper. 



A. Age problem 

Some cosmological models suffer from the so-called "age problem" , namely, the existence 
of high redshifts objects whose age at some redshift seem to exceed the Universe's age 



predicted at that redshift (as in the ACDM model, see e.g. |24|). 
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The age of the Universe at redshift z is 



t{z) = t 



dz' 



(l + z')H(z>) 



(21) 



Figure H] shows the age of the Universe as a function of redshift for the best fit holographic 
model and the ACDM. Also marked in the figure are the ages and redshifts of three luminous 
old objects: galaxies LBDS 53W069 (z = 1.43, t = 4.0 Gyr) |25| and LBDS 53W091 



{z = 
Gyr) 



.55,1= 3.5 Gyr) 



26 



24 



27J, as well as the quasar APM 08279+5255 (z = 3.91, t = 2.1 



28j. As is apparent, the ages of two first old objects result compatible with both, 



the holographic and the ACDM model; however, the age of the old quasar falls only within 
2cr with the ages predicted by these two models. Thus, some tension exists in this regard. By 



contrast, the interacting holographic model of Ref. [23(, which takes as infrared cutoff the 
Hubble radius, is free of the problem. At any rate, it remains to be seen in which direction 
(if any) future observations will "move" the age of the said quasar. 



M 
>i 





Figure 4. Age of the Universe, t, in Giga-years vs. redshift. The solid line corresponds to the best 
fit holographic model and the dashed line to the ACDM model. The data points from left to right 
locate the old objects LBDS 53W069, LBDS 53W091 and APM 08279+5255. 
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III. DISCUSSION OF THE COSMIC COINCIDENCE 



In holographic models of dark energy that take the Ricci's length as the infrared cutoff 
one can obtain a finite and approximate constant ratio r for an ample redshift span even if 
no interaction between the dark components is assumed -see e.g. |7j. Here we analyze how 
this comes about, and we note that while this approach alleviates the coincidence problem 
it does so only partially since it would entail that we are living in a special time. 

We start by rewriting Friedmann's equation (J2J) with the help of the saturated holographic 



bound, p x = 3M 2 c 2 (H + 2H 



as 



3H 2 = 3c 2 {H + 2H 2 ) + 



^(1 + z) 3 
M 2 [ +Z) ' 



(22) 



and, for convenience, introduce the ancillary variable y 
form 



3a 



2 y_ 

y 2 





y 








jj 


y 2 _ 





+ 



1 + z. Thus, Eq. fl22J takes the 

(23) 



Pmo 
M 2 



y 



-3a 



By equating coefficients, we get a = c 2 /(2c 2 — 1), and 



y + 



Pmo y-3a+i 



3c 2 M 2 a 



0. 



(24) 



Multiplying the latter by y the differential equation can be readily solved. Upon reverting 
to the original variable we obtain 

2 



3H Z 



(2 - c 2 )M\ 



-p M0 (l + zY + (3Mp 2 p M0 (l + z) 2 ^ 



(25) 



where /3 is a positive-definite integration constant that can be identified as /3 
and, of the order of unity since ro and c 2 lie not far from 0.4. 



2-c 5 



Recalling Eq. (0) we finally get 



Px 



r p M0 (l + z? + (3 Pmo (I + zf 



(26) 



2-c 2 ' 

The DE density is contributed by two terms. The first one redshifts exactly as non-relativistic 
matter. The second one, in view that c 2 is bounded by 0.36 < c 2 < 0.8, results subdominant 
for z of order of unity and larger. Therefore, we can safely conclude that r = Pm/px — (2 — 
c 2 )/c 2 for < z. This is why the ratio r results of order unity in an ample redshift interval, 



also in the absence of interaction, as in Ref. 



7|. However, in view of the observational lower 
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bound on c 2 , we see that r — > as z — > — 1. So, in the holographic non-interacting model, 
r results well below unity, close to zero, and approaches this null value asymptotically for 
an infinite span of time. Altogether, according to this model, we live in a very special and 
transient period in which r results comparable to unity. 

IV. OBSERVATIONAL CONSTRAINTS 

To constrain the four free parameters (flxo, c 2 , rf, and H ) of the holographic model 
presented above we use observational data from SN la (557 data points), the CMB-shift, 
BAO, and gas mass fractions in galaxy clusters as inferred from x-ray data (42 data points), 
the Hubble rate (15 data points), and the growth function (5 data points). Being the 
likelihood function defined as C oc exp(— x 2 /2) the best fit follows from minimizing the sum 

Xtotal Xsn Xcmb "X-bao "X-x—rays ^Hubble Xgf- 



A. SN la 

We contrast the theoretical distance modulus 

fi th (zi) = 51og 10 fy^J + tM) , (27) 
where /x = 42.38 — 51og 10 h, with the observed distance modulus ii bs{ z i) of the 557 SN la 



compiled in the Union2 set [29|. The latter assemble is much richer than previous SN la 
compilations and has some other advantages, especially the refitting of all light curves with 
the SALT2 fitter and an enhanced control of systematic errors. In (12"7|) D L = (1+z) Jq E(z'-p) 
denotes the Hubble- free luminosity distance, with p the model parameters (flxo, c 2 , r/, and 
H ), and£(z;p) := H{z;p)/H . 

The x 2 from the 557 SN la is given by 



2 ,' \ \P>th(Zi) — Vobs(Zi)] 2 / oq n 



i=l v 



where cij denotes the la uncertainty associated to the ith data point. 



To eliminate the effect of the nuisance parameter /i we resort to the method of 30j to 
obtain xL = x 2 st inimum) = 543.70. 
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B. CMB-shift 



The displacement of the first acoustic peak of the CMB temperature spectrum with 



respect to the location it would take shou 
Sitter model is given by the CMB-shift 



31 



d t 



re Universe be described by the Einstein-de 



32| 



K 



dz 



(29) 



where z rec ~ 



1089 is the redshift at the recombination epoch. This parameter is approx- 
imately model-independent but not quite as the above expression somehow assumes the 
ACDM model. 

The 7-year WMAP data provides TZ(z rec ) = 1.725 ± 0.018 The best fit value of the 
model is TZ(z rec ) = 1.727 ± 0.030. Minimization of 



XcmbVP) 



(TZ 



th 



TZ 



obs I 



a 



(30) 



TZ 



i 2 (mirvi't 

produces Xcmb- 



2 (minimum) 
shift 



0.013. 



BAO 



Pressure waves originated from cosmological perturbations in the primeval baryon-photon 
plasma produced acoustic oscillations in the baryonic fluid. These oscillations have been 
unveiled by a clear peak in the large scale correlation function measured from the luminous 
red galaxies sample of the Sloan Digital Sky Survey (SDSS) at z — 0.35 |34| as well as in the 
Two Degree Field Galaxy Redshift Survey (2dFGRS) at z = 0.2 [35]. These peaks can be 
traced to expanding spherical waves of baryonic perturbations with a characteristic distance 
scale 

ZBAO_^j Z J A ° dz 



DvizBAo) 



-see e.g. 



36] 



H(z BAO ) 



H(z) 



(31) 



Data from SDSS and 2dFGRS observations yield D v (0.35)/D v (0.2) = 1.736 ±0.065 |35|. 
The best fit value for the holographic model is -D„(0.35)/Z\(0.2) = 1.664 ± 0.003, and 
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minimization of 



2 ( v _ ([£)„(0.35)/D„(0.2)] f/t - [^(0.35)/^(0.2)] ofes ) 2 

a D v (0.35)/D v (0.2) 



2 (minimum) , 

§ lves Xiao = 1-20. 



D. Gas mass fraction 

As is well known, a very useful indicator of the overall cosmic ratio ^ibaryons/^M, nearly 
independent of redshift, is the fraction of baryons in galaxy clusters [37] ■ This quantity can 
be determined from the x-ray flux originated in hot clouds of baryons and it is related to the 
cosmological parameters by f gas oc d 3 / 2 , where d A '■= (1 + z)~ l J* denotes the angular 
diameter distance to the cluster. 

We used measurements by the Chandra satellite of 42 dynamically relaxed galaxy clusters 
in the redshift interval 0.05 < z < 1.1 [38]. In fitting the data we resorted to the empirical 
formula 

fM ^^^(<Pf (33 , 
1 + s(z) il M0 \ d A J 



(see Eq. (3) in Ref. 38]) in which the ACDM model served as a reference. Here, the 
parameters K, A, 7, b(z) and s(z) model the aboundance of gas in the clusters. We set 
these parameters to their respective best fit values of Ref. \&\ . 
The x 2 function from the 42 galaxy clusters reads 

E ([fgas( z i)]th ~ [J gasK z i)\obs) /„ ,\ 

/T-2 ( V , I 



2 i 1 \ ^ {[fgas{ z i)]th [fgas{ z i)]obs) 

/(.x— raysvP) . . , . 



1 • ; • 1 ; 1 2 (minimum) , -, 

and its minimum value comes to be Xx-rays — 41.79. 
Figure [5] shows the fit to the data. 



E. History of the Hubble parameter 



Recently, high precision measurements by Riess et al. at z = 0, from the observation of 
240 Cepheid variables of rather similar periods and metallicities 39| , as well as measurements 
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Figure 5. Gas mass fraction in 42 relaxed galaxy clusters vs. redshift. The solid and dashed curves 
correspond to the best fit models: holographic and ACDM, respectively. The data points with 
their error bars are taken from Table III of Ref . 38] . 



by Gaztanaga et al, at z = 0.24, 0.34,and0.43 [40j, who used the BAO peak position as a 
standard ruler in the radial direction, have somewhat improved our knowledge about H(z). 
However, at redshifts above, say, 0.5 this function remains largely unconstrained. Yet, in 
order to constrain the holographic model we have employed these four data alongside 11 
noisier data in the redshift interval 0.1 ^ z ^ 1.8, from Simon et al. 4l| and Stern et al. 



42J, obtained from the differential ages of passive-evolving galaxies and archival data. 
By minimizing 



XHubbleVP) 



15 

£ 

1=1 



[Hthjzj) - H obs { Zi )f 
ctHzA 



(35) 



we got XnubWe" 1 ""^ = 9-^7 and H = 71.8 ± 2.9 km/s/Mpc as the best fit for the Hubble's 
constant. Figure |6] depicts the Hubble history according to the best fit holographic model 
and the best ACDM model. 



F. Growth function 

Up to now, we have considered background quantities that chiefly depend on H(z) whence 
they are not very useful at discriminating between cosmological models that present a similar 



14 



25C 



200 



15C 



e 100 
2 



50 



£2 X0 =0.707 r f =0.013 c =0.407 



la 

la Ac dm 

AC DM 




-0 . 5 



Figure 6. History of the Hubble factor in terms of the redshift for the best fit values of the 



holographic model (solid 
are borrowed from Refs. 



ine) and the ACDM model (dashed line). The data points and error bars 



39|442] 



Hubble history, independently of how different they are otherwise. One manner to circum- 
vent this hurdle is to study evolution of the growth function / = din Sm/ din a, where Sm 
denotes the density contrast of non-relativistic matter. 

The evolution of the latter obeys the coupled set of equations 

Sm ~ — vm = -— —QHS M , (36) 
a r 

v M + Hv M + -0 = 0, (37) 
a 

where the Newtonian potential fulfills Poisson's equation 

k 2 

— = -AttGpmSm . (38) 
cr 

Solving the equations, and expressing in terms of r and using w = [c 2 (l + r) — 2]/(3c 2 ) 
obtained from PJ, c _2 (l + r) _1 = (H / H 2 ) + 2, and (TTTj) . we get from the evolution equations 
(ED and ([10]) 



w , f 2 , ( 1 , . 3r 3 + 2g 2 (l + r) 3 -2gr(2 + r-r 2 ) _ n 
; +/ + Q ^~) f 2r2(l + r) -°' (39) 

where /' = df/dlna. Note that in the limit Q — > 0, last equation collapses to the corre- 
sponding expression of the Einstein-de Sitter scenario (Q M (z) = 1 and Sm oc a oc i 2 / 3 ); that 
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is to say, f' + f 2 + [2 + (H/H 2 )}f = 30 M /2. (Recall that in Einstein-de Sitter H/H 2 = -3/2 
and the solution of the equation for / is simply / = 1). In constraining the model we have 
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Figure 7. Growth function, /, vs. redsfiift as follows from integrating Eq. (|39p (solid line). Also 
shown is the prediction of the ACDM model (dashed line). The observational data were taken from 
Ref. [43|]. In constraining both models we have used only the five less noisy data depicted in the 
figure, (i.e., data corresponding to redshifts below 1.5). 



taken only the five lowest redshift data of the growth function shown in Fig. [7] -the other 
data present very large error bars. The best fit yields y 2 = 1.06. 



Figures IBl - fTUl and table I summarize our findings. 

Figure [TOl depicts the la confidence contours from SN la (dashed yellow), CMB-shift (solid 
black), BAO (dashed blue), x-rays (dashed black), history of the Hubble function (dot-dot 
dashed green), and grow function (dot-dashed red) in the (Qxo,c 2 ) plane (left panel) and 
the (Q X o, r f) plane. The joined constraints corresponding to Xtotal are shown as shaded 
contours. As is apparent from left panel most of the discriminatory power arises from the 
near orthogonality between the x-ray and CMB-shift and supernovae contours. However, in 
the right panel the supernovae contour appears nearly degenerated with respect to the x-ray 
contour. 
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Figure 8. Panels from left to right show the 68.3% and 95.4% confidence contours for the pair 
of free parameters (fixo? Hq), (Qxo, c 2 ), (Qxq, rf), respectively, obtained by constraining the 
holographic model with SN Ia+CMB-shift+ BAO+x-ray+H(z)+growth function data. The solid 
point in each panel locates the best fit values. 
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Figure 9. Same as Fig. [8] for the pairs of free parameters (Hq, rf), (c 2 , r/), and (Hq, c 2 



Altogether, by constraining the holographic model of Section II with data from SN la, 
CMB-shif, BAO, x-rays, H(z), and the growth function we obtain Qxo — 0.707 ± 0.009, 
c 2 = 0.407±g;g||, r f = 0.013+°;°^, and H = 71.8 ± 2.9 km/s/Mpc as best fit parameters. It 
is worth noticing that the non- interacting case, Q = (which implies rj = via Eq. ([12]) ). 
lies over 2cr away from the best fit value. This feature seems typical of holographic dark 



energy models (see e.g. 



7. 9. 



23|) 



Table [T] shows the partial, total, and total y 2 over the number of degrees of freedom of 
the holographic model along with the corresponding values for the ACDM model. In the 
latter one has just two free parameters, Qmo an d H . Their best fit values after constraining 
the model to the data are fi A /o = 0.266 ± 0.006, and H = 71.8 ± 1.9km/s/Mpc. 

Although the ACDM model fits the data somewhat better, A% 2 ~ 2.5, than the holo- 
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Figure 10. Left panel: The probability contours for SN la, CMB, BAO, x-ray, H(z), and growth 
function, in the (fl xo , c 2 ) plane. The joined constraint uses x? otal = Xsn + x\ m b + xL> + xl-rays + 
XHubbie + X 2 g f Right panel: Idem for the (n xo ,r f ) plane. 



Model 


Xsn 


Xcmb 


Xbao 


Xx— rays 


Xh 


4f 


Xtotal 


xlotJdof 


Holographic 


543.70 


0.01 


1.20 


41.79 


9.57 


1.06 


597.34 


0.96 


ACDM 


542.87 


0.05 


1.13 


41.59 


8.73 


0.43 


594.80 


0.96 



Table I. x 2 values of the best fit holographic model (tt x = 0.707 ± 0.009, c 2 = 0.407±g;g||, r f = 
0.013j#ooi, and H = 71.8 ± 2.9 km/s/Mpc), and the best fit ACDM model (O M0 = 0.266 ±0.006, 
and Hq = 71.8 ± 1.9 km/s/Mpc). 



graphic model, that has two more free parameters, it its uncertain which model should be 
preferred in view that the former cannot address the cosmic coincidence problem and the 
latter substantially alleviates it. More abundant and accurate data, especially at redshifts 
between the supernovae range and the CMB, will help decide the issue. Nonetheless, we 
believe that the uncertainty will likely persists till a breakthrough on the theoretical side 
allows us to calculate with confidence the true value of the cosmological constant. 



V. CONCLUDING REMARKS 



We performed a statistical study of the best fit parameters of the holographic model 
-presented in section [TTJ- using data from SN la, CMB-shift, BAO, x-ray, the Hubble history, 
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and the growth function; 621 data in total. The maximum likelihood (or minimum \ 2 ) 
parameters are tt xo = 0.707 ± 0.009, c 2 = 0.407±g;g||, r f = 0.013±g;gg|, and H = 71.8 ± 2.9 
km/s/Mpc with x 2 /dof ~ 0.96. The fi™ and H Q values fall within la of the corresponding 
values determined by Komatsu et al. 33j (0.734 ± 0.029 and 71.0 ± 2.5 km/s/Mpc, respec- 
tively). The evolution of the equation of state parameter, w, at redshift below 1.2 (see right 
panel of Fig. [2]) is compatible with the observational constraints derived in 21[ and, as in 
other Ricci's holographic models 7H10|. it crosses the phantom divide line at recent times 
(see right panel of Fig. [2]). Curiously enough, the flxo and c 2 best fit values of this model 
agree within la with the corresponding values obtained by Suwa et al. [9j despite the use 
of a very different interaction term between the dark components. 

It is to be emphasized that holographic models do not contain the ACDM model as a 
limiting case since the energy density of the quantum vacuum, being constant, cannot be 
holographic. It is also noteworthy that, in general, the c 2 term in the holographic expression 
for the dark energy, Eq. (CQ), should not be considered constant except precisely when the 
Ricci's length is chosen as the infrared cutoff 44]. 

A lingering problem, both for this model and the ACDM model, refers to the age of the 
old quasar APM 08279+5255 at redshift z = 3.91. In both models the measured quasar 
age would fall within la only if the Hubble constant, Hq, would come down substantially 
-something we do not expect though it cannot be excluded. Accordingly, we must wait for 
further observational data to see whether the present tension gets exacerbated or disappears. 

Before closing, it results interesting to contrast the model explored in this paper with 
the model of Ref. [23)]. The present one shows a better fit to the growth function at low 
redshifts, as well as to the CMB-shift and BAO data. However, it does not fit so well the 
age of the old quasar APM 08279+5255. 
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